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Abstract.  Thermophoresis  of  axially  symmetric  bodies  is  investigated  to  first  order  in  the  Knudsen-number, 
Kn.  The  study  is  made  in  the  limit  where  the  typical  length  of  the  immersed  body  is  small  compared  to  the 
mean  free  path.  It  is  shown  that  in  this  case,  in  contrast  to  what  is  the  case  for  spherical  bodies,  the  arising 
thermal  force  on  the  body  is  not  in  general  antiparallel  to  the  temperature  gradient.  It  is  also  shown  that 
the  gas  excerts  a  torque  on  the  body,  which  in  magnitude  and  direction  depends  on  the  body  geometry. 
Equations  of  motion  describing  the  body  movement  are  derived.  Asymptotic  solutions  are  studied. 


INTRODUCTION 

Thermophoresis  as  a  phenomenon  has  been  known  for  a  long  time,  and  several  authors  have  approached 
the  problem.  For  example,  Einstein  calculated  the  final  velocity  of  a  spherical  particle  in  a  heat  conducting 
gas  using  elementary  kinetic  theory.  A  recent  review  of  the  phenomenon  is  given  in  an  article  by  Sone,  see 
[1].  Another  review  is  given  by  Talbot,  see  [2].  The  first  systematic  attempt  to  describe  the  thermophoresis 
phenomenon  using  kinetic  theory  is  found  in  an  article  by  Waldmann  from  1959,  see  [4] .  That  work  was  made 
under  the  assumption  that  the  mean  free  path  of  the  gas  is  much  larger  than  the  dimension  of  the  body.  Further 
results  are  found  in  a  variety  of  articles,  but  these  mostly  apply  to  larger  bodies,  using  asymptotic  methods, 
and  often  to  spheres,  cf  [3]. 

In  this  work,  thermophoresis  of  axially  symmetric  bodies  is  considered.  The  typical  macroscopical  length  L 
over  which  the  temperature  varies  is  assumed  to  be  much  larger  than  A,  the  mean  free  path  of  the  gas.  To 
this  end  we  define  the  Knudsen  number  according  to  Kn  =  X/L.  Thus  Kn  <C  1,  and  we  can  use  the  first  order 
Chapman-Enskog  solution  to  the  Boltzmann  equation,  see  [5],  to  describe  the  heat  conducting  gas  from  the 
macroscopical  viewpoint.  The  typical  diameter  R  of  the  body  is  assumed  to  bee  much  smaller  than  A.  R  is 
also  assumed  to  be  much  larger  than  d,  the  diameter  of  a  gas  molecule.  In  calculating  the  momentum  transfer 
from  the  gas  to  the  body  we  therefore  use  the  equations  of  free  molecular  flow  (  [5]).  We  also  assume  that 
Ma  =  v/cs  <C  1,  where  v  is  the  speed  of  a  body  surface  element  and  where  cs  is  the  speed  of  sound. 


LOCAL  MOMENTUM  TRANSFER  TO  THE  BODY 

If  the  one-particle  distribution  function  /  =  /(x,  c,  t )  describing  the  gas  is  known  and  the  proper  gas-body 
surface  boundary  conditions  are  specified  it  is  in  principle  possible  to  calculate  the  net  transfer  of  momentum 
from  the  gas  molecules  to  a  body  immersed  in  the  gas.  Let  us  first  of  all  consider  the  interaction  between 
the  gas  and  a  surface  element  of  the  body,  with  area  dS  and  normal  n,  that  is  at  rest  in  the  gas.  The  net 
momentum  transfer  from  the  gas  to  the  body  surface  element  can  now  be  expressed  by 

dFi  =  —  f  f  mciC  ■  n f^d3c+  f  mciC  ■  nf^d3c  \  dS  (1) 

c*n<0  J  cn>0  / 

Here  /W  describes  the  incident  molecular  stream  and  f(r}  describes  the  reflected  stream,  m  is  the  mass  of  a 
gas  molecule,  and  c  is  the  velocity.  The  surface  element  is  taken  to  be  impermeable,  and  thus  the  stream  of 
gas  molecules  incident  to  the  surface  element  equals  the  stream  of  outgoing  molecules.  Therefore 
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c  •  n f^d3c  =  0. 


(2) 


[  cn /«d3c  +  [ 

J  c*n<0  J  cn>0 

As  the  typical  length  of  the  body  is  small  compared  to  the  mean  free  path  in  the  gas  and  as  the  body  surface 
is  taken  to  be  convex,  the  incident  stream  of  molecules  is  unaffected  by  the  presence  of  the  body,  and  thus 
/(i)  _  j  where  /  is  the  distribution  function  for  the  gas  in  the  absence  of  the  body.  The  reflected  stream 
is  determined  from  a  Maxwell-type  boundary  condition  and  aquires  thus  upon  reflection  two  separate  parts: 
One  part  that  is  specularly  reflected,  and  one  part  that  is  diffusely  reflected.  Here,  the  relative  amount  of 
momentum  tangential  to  the  surface  element  transferred  to  the  body  from  the  incident  stream  is  measured  by 
the  momentum  accommodation  coefficient  aT,  see  [5].  This  boundary  condition  has  the  following  form: 

f(r)  =  (1  -  ar)f  +  aT(n«/n(°))/f°),  (3) 

where  f  =  /(c'),  in  which  c'  =  c,;  —  ‘IniiijCj.  is  thus  the  distribution  function  of  a  specularly  reflected  stream. 
The  second  term  describes  the  diffusely  reflected  part  of  the  outgoing  stream.  is  a  Maxwellian  distribution 
function  given  by 


f(0) 


=  n 


(27rkBTys/2  c 


(4) 


describing  gas  molecules  which  are  in  thermal  equilibrium  with  the  surface  element.  In  this  expression,  nr  is 
the  number  density  of  the  reflected  stream,  and  the  dimensionless  velocity  C  is  defined  by 


V  m 

T  is  the  local  gas  temperature  and  kB  is  the  Boltzmann  constant.  It  is  assumed  here  that  the  temperature  of 
the  body  is  the  same  as  the  local  temperature  of  the  gas.  The  integral  (1)  can  now  be  simplified.  Introducing 
a  dimensionless  mass  flux  and  a  dimensionless  momentum  flux  incident  on  the  body  according  to 


M.  =  — —  f  CiUifSC  and  Mf  =  --  '  f  CiCjrijftfC 

n  \  m  J  JCini<0  n\  m  J  JCini<o 


(5) 


the  force  on  a  surface  element  becomes,  where  p  =  nkBT, 


dFi  =  p  [4(1  —  aT)nirijMj  +  aT  (2Mi  —  sfn Mrii')]  dS.  (6) 

In  the  present  case,  we  shall  consider  a  body  which  is  small  compared  to  the  mean  free  path.  This  body 
will  start  to  move  under  the  influence  of  the  force,  and  thus  the  surface  element  will  experience  a  homogenous 
flow  of  the  surrounding  gas.  This  flow  will  excert  an  additional  force  on  the  surface  element.  Therefore,  when 
calculating  the  force  on  the  surface  element  according  to  (6),  we  consider  a  distribution  function  corresponding 
to  a  gas  which  is  subject  to  a  temperature  gradient  and  a  homogenous  flow.  It  is  assumed  that  Kn  =  A/A, 
where  A  is  the  mean  free  path  and  where  L  is  the  macroscopic  length  scale,  is  small.  We  also  assume  that  the 
resulting  motion  of  the  surface  element  is  small  compared  to  the  speed  of  sound. 

For  a  gas  subject  to  a  non-uniform  temperature  distribution  the  distribution  function  takes  the  form,  see  [5], 


/  =  /o  (1  +  4>) , 


(7) 


where  (f>  can  be  found  through  the  Chapman-Enskog  expansion  to  first  order  in  Kn.  The  result  for  pure  heat 
conduction  is 


The  function  A  is  usually  expressed  in  Sonine  polynomials  according  to 


(8) 


A(C2)  =  —  ^2  anS^l(C2),  with  m 

n—1 


2  mn 

5k%T' 


(9) 


868 


Here  n  is  the  heat  conductivity.  If  we  introduce 

A(C2)  =  A{C2)/au 


we  may  express  <fi  in  terms  of  the  heat  current  q,  =  — kT*  according  to 

,  4  / 2kBT\ _1/2 

A(c 

Now  <f>  is  integrated  to  give  the  moments  (5)  and  one  arrives  at 


ka  1  a  k  a  1  ,  7  f2wkBT\ 

M  =  and  Mi  =  --Tit  +  —  - 

2y/-K  4  lOp  \  m  J 


-1/2 


{Sij  +  riiTij)  qj. 

The  formulas  (11)  thus  apply  to  a  body  surface  element  at  rest.  Here,  the  integral  7  is  given  by 


7  = 


pOO 

/  C5A(C2)e 
Jo 


dC. 


(10) 


(11) 


(12) 


If  only  the  first  term  in  the  Sonine  polynomial  expansion  of  A  is  retained,  7=1.  The  corresponding  integral 
in  M, 


dC, 


vanishes  since  the  gas  is  at  rest. 

The  temperature  gradient  will  give  rise  to  a  force  on  the  immersed  body  and  thus  set  the  body  in  motion. 
Following  Waldmann,  if  we  denote  the  velocity  of  a  body  surface  element  by  u-n  this  motion  is  taken  into 
account  by  transforming  the  distribution  function  to  a  frame  of  reference  in  which  the  surface  element  of  the 
body  is  momentarily  at  rest.  When  the  transformation  Cj  — >  c*  —  Ui  is  made,  both  the  Maxwellian  part  and 
the  heat  conduction  part  of  (7)  is  transformed.  However,  since  Ma  is  small,  the  cross  term  coupling  translation 
and  heat  conduction  will  be  of  the  order  of  KnMa,  and  will  thus  be  ignored  in  the  present  context.  The 
moments  M  and  Mi  get  additional  contributions  of  the  order  of  Ma,  and  the  result  is 


M  = 


(2irkBTy1/2 


TIi'U'i  5 


(13) 


Mi  =  \  (27r„f 7 )  + ninJ )  •  (14) 

The  resulting  force  on  a  body  surface  element  (6)  becomes 

dFi  =  -priidS  +  (^~^j  { ar  +puij  +  Q(4  -  3 aT)qj  +p  [4  -  (3  ~  «r]  rijn^dS  (15) 

This  force  will  in  general  also  produce  a  torque  acting  on  the  body  through  the  surface  element  according  to 

—  f  ij  /,-  x j  d  /'/- ,  ( 1  ii ) 

where  x  is  the  vector  from  the  center  of  mass  of  the  body  to  the  surface  element.  As  the  force  and  the  torque 
are  applied  to  the  body,  a  velocity  and  also  an  angular  velocity  will  be  produced.  Therefore,  the  appropriate 
Gallileian  transformation  to  be  used  in  this  context  is 

Ui  =  (Vi  -f-  CijkLdjXk'), 

where  the  minus  sign  stems  from  transforming  to  the  frame  of  reference  where  the  gas  is  at  rest.  The  integration 
of  the  force  and  the  torque  over  the  total  body  surface  involves  calculating  geometrical  tensor  integrals  of  the 
type 

(17) 

These  are  all  isotropic  functions  of  the  unit  vector  N  defining  the  axis  of  symmetry,  and  are  listed  in  the 
appendix. 


L 


nilni2...nikxhxh...xjldS  =  A 
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FORCE  AND  TORQUE  ON  THE  BODY 


Now  we  are  in  a  position  to  write  down  the  expressions  for  the  force  and  the  torque  acting  on  the  body.  In 
what  follows,  the  dimensionless  dynamical  coefficients  ai,..,ag  are  introduced  for  simplicity.  They  consist  of 
products  between  linear  functions  of  aT  and  geometrical  coefficients  given  by  contractions  of  the  geometrical 
tensor  integrals  encountered  in  the  previous  section(17).  The  coefficients  a\, ...  ag  are  listed  in  the  appendix. 

The  force  becomes 

Fi  =  |  (sij  +  aiN<iNj> )  Qj  -  Ps  («2%  +  a3iV<iiVJ>)  Vj  +  3pS3/2a4eijkNjujk^  .  (18) 

The  first  term  on  the  right  hand  side  of  (18)  is  a  thermal  contribution  to  the  force  proportional  to  the  heat 
current.  It  contains  one  isotropic  part,  a.^Sij  and  one  symmetric  traceless  part,  (I2N ^jN 3> .  While  the  former 
part  acts  on  a  spherically  symmetric  body,  the  latter  one,  which  vanishes  in  case  of  spherical  body  surface, 
projects  momentum  onto  the  direction  of  the  axis  of  symmetry.  The  order  of  magnitude  of  this  term  is  ~  pSKn. 
The  second  contribution  to  the  total  force  is  a  friction  force  with  geometrical  features  similar  to  the  thermal 
force.  The  last  term  is  a  force  acting  on  the  body  as  a  consequence  of  the  body  rotation.  The  two  last  terms 
both  have  the  magnitude  ~  pSMa. 

The  torque  is  given  by 

^  U5 ^ijkNjQk  3 pS  ^  Q,Q€ijk,N jVfc  “  ( ^7  A/  T  Ug / P\j y. )  UJj  ^  .  (19) 

In  the  expression  for  the  torque,  the  first  term  is  a  thermal  torque  that  acts  to  align  the  axis  of  symmetry  to 
the  heat  current.  This  contribution  is  of  the  order  of  ~  pS3^2Kn.  The  second  term  is  active  whenever  the  axis 
of  symmetry  is  not  parallel  or  antiparallel  to  the  velocity.  The  third  term  is  a  friction  torque.  The  two  latter 
contributions  are  of  the  order  of  VS3/2Ma. 

EQUATIONS  OF  MOTION 

In  order  to  formulate  the  equations  of  motion  for  the  body  immersed  in  the  heat  conducting  gas,  we  introduce 
a  body-fixed  orthonormal  frame  of  reference  with  the  origin  in  the  body’s  center  of  mass.  The  basis  vectors  are 
denoted  by  e(a)  where  the  index  a  =  1,2,3  numbers  the  basis  vectors.  One  of  these  basis  vectors  is  naturally 
chosen  to  coincide  with  the  axis  of  symmetry  N.  The  equations  of  motion  become 

)§ ( mBVi )  =  (277^fr)  1/2  |  ( 6ij  +  aiIV<ilVJ>)  qj  -  pS  ( a28ij  +  agN^Njy)  vj  +  3pS3/2aieijkNjujk  j- , 

ftlhjUj)  =  )-x/ 2  ^a5£ijkNjqk  _  3pS3/2  a6eijkNjVk  -  (a7^  +  asN^N^)  Uj}  , 

fMa) )  =  eb'fcwi4a)> 

(20) 


where  is  the  body  mass  and  where  is  the  inertia  tensor.  This  tensor  has  the  form 


hj  —  To  ( biSij  +  b2N<iNj>) . 


(21) 


For  values  of  see  Appendix. 

It  is  possible  to  construct  two  asymptotic  solutions  to  these  equations  by  observing  that  F  =  M  =  0  if  v 
is  parallel  to  q  and  N  is  parallel  or  antiparallel  to  q,  together  with  u;  =  0:  The  equation  of  motion  is  then 
reduced  to  an  expression  for  the  asymptotic  velocity  V: 


4  3  +  2ai  _  7  1  +  \  (4  -  3ar)  c^2,0) 

15^3ci2  +  2a3  5p  1  +  +  \  [8  —  (6  —  7r)  aT]  c(2>°)  ^ 


(22) 


Here,  7  is  given  by  (12).  It  depends  on  the  intramolecular  forces  of  the  gas  molecules;  7  =  1  if  only  the  first 
term  of  the  Sonine  polynomial  expansion  of  the  heat  conducting  part  of  the  distributionfunction  is  retained. 
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The  new  terms  are  those  containing  the  parameter  c(2,0)  that  depends  on  the  body  geometry.  c(2,°)  measures 
deformations  of  the  body  surface  from  a  sphere:  For  a  sphere,  c^2jr>  =  0.  In  this  case  the  asymptotic  velocity 
(22)  coincides  with  Waldmann’s  result,  [4].  If  the  the  body  is  extremely  oblate,  or  shaped  like  a  coin,  c*2,0)  =  1. 
If  the  body  is  extremely  prolate,  or  shaped  like  needle,  r/2,0)  =  —  -The  order  of  magnitude  of  this  velocity 
is  Kn  ■  cs. 

Non-dimensional  variables  are  now  introduced  according  to  V{  =  t^v*,  t  =  ^  (27rfrf T)  t*,  and 

(2,r^fr)  w*,  where  q=  |q|.  The  non-dimensional  heat  current  is  chosen  to  be  q  =  (1,0,0).  A 
simulation  of  the  non-dimensional  equations  of  motion  yielding  the  body  orbit  of  a  cone  and  a  coin  is  plotted 
in  figure  1. 
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FIGURE  1.  Orbit  of  a  cone  and  a  coin,  with  the  heat  current  in  the  ^-direction.  Initially  the  axis  of  symmetry  N 
of  the  two  bodies  makes  the  angle  45°  with  the  heat  current.  The  initial  non-dimensional  velocity  of  the  two  bodies  is 
5  and  is  directed  along  N.  The  spacing  between  the  points  representing  the  body  orbits  corresponds  to  a  fixed  time 
interval.  Note  that  the  cone  assumes  a  final  velocity  parallel  to  the  heat  current,  whereas  the  final  velocity  of  the  coin 
also  has  a  component  in  the  plane  orthogonal  to  the  heat  current. 


Numerical  simulations  seem  to  suggest  that  the  bodies  reach  asymptotic  states  after  a  distance  of  the  order 
of  1  in  non-dimensional  units,  that  is,  l  =  (27r^f T)  ~  T ^  dimensional  units.  Here,  R  is  the 

dimension  of  the  body,  and  d  is  the  diameter  of  a  gas  molecule.  The  difference  in  asymptotic  velocity  is  due  to 
a  symmetry  in  the  body  shape  of  the  coin  that  is  lacking  in  the  body  shape  of  the  cone.  This  is  explained  in 
the  proceeding  paragraph. 

In  order  to  investigate  the  stability  of  the  two  asymptotic  states  we  linearize  the  equations  of  motion  around 
the  state  where  N  is  nearly  parallel  to  q.  The  eigenvalues  of  the  resulting  linear  system  is  calculated  for  a 
special  type  of  body,  a  ’double  cone’,  shown  in  figure  2.  This  body  consists  of  two  cones  with  a  common  base. 
The  radius  of  the  base  is  denoted  by  D,  and  the  total  length  by  R.  The  base  is  situated  a  distance  s  •  R  from 
the  left  cusp,  where  the  dimensionless  parameter  obeys  0  <  s  <  1.  When  s  =  0  the  double  cone  degenerates  to 
a  single  cone  with  its  cusp  pointing  in  the  direction  of  N.  When  s  =  1  we  recover  another  single  cone,  pointing 
in  the  direction  of  — N. 


R 


Plotting  the  eigenvalues  against  s  yields  a  picture  of  the  preferred  orientation  of  the  body.  It  turns  out  that 
stabillity  is  determined  by  two  identical  eigenvalues  (the  other  eigenvalues  are  either  negative  or  zero  for  all 
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values  of  s).  This  eigenvalue  is  plotted  against  s  in  figure  3  for  a  somewhat  prolate  body  ( R/D  =  3)  and  for 
an  oblate  body  ( R/D  =  1/3). 


FIGURE  3.  Critical  eigenvalue  plotted  against  s  for  two  different  double  cones 

It  is  apparent  that  the  stable  state  for  the  prolate  body  occures  when  s  <  0.5  and  for  the  oblate  body  when 
s  >  0.5.  This  may  indicate  that  the  stabillity  is  determined  by  the  orientation  of  the  center  of  mass  of  the 
body  and  the  center  of  mass  of  the  surface  of  the  body  along  the  direction  of  the  heat  current. 


BODY  INVARIANT  UNDER  REFLEXION  IN  A  PLANE  ORTHOGONAL  TO 

THE  AXIS  OF  SYMMETRY 


Now  we  consider  a  body  that  is  symmetric  with  respect  to  an  equator,  that  is,  a  body  that  is  invariant 
under  the  transformation  N  — >  — N.  Then  the  coefficients  0,4, 125  and  ae  vanish.  It  it  easy  to  see  that  then  the 
rotational  motion  decouples  from  the  translation  and  becomes  independent  of  the  heat  current.  The  remaining 
part  of  the  angular  momentum  equation  in  (20)  is  then 

d  fr  ,  pS2  f2nkBT\~1/2  .  c  .  /00, 

^(IijUj)  —  2  (  m  )  (orSij+atNciNjJuj.  (23) 

This  equation  states  that  as  t  — >  00  all  constant  values  of  N  is  a  solution.  In  the  momentum  equation  (20)  the 
coupling  to  u)  vanishes,  but  the  dependence  on  N  still  remains  through  the  tensor  N<iNj>.  The  asymptotic 
solution  to  the  momentum  equation  in  this  case  becomes 

Vi  =  [Sij+k(6ij-NiNj)}Vj,  (24) 

where  Vj  is  the  asymptotic  velocity  (22)  found  in  the  previous  section.  The  tensor  Sij  —  N{  Nj  picks  up  the 
projection  of  q  that  lies  in  the  plane  perpendicular  to  N.  The  parameter  k  is  introduced  for  simplicity  and  is 
given  by  (we  recall  that  r/2,0)  is  a  geometrical  parameter  encountered  in  (22)) 


2(8  +  ttov)  +  [16  -  (12  -  27r)ar  -  37ra2]c(2,o)  _  [32  _  (48  _  47r)aT  +  (is  -  3tt)q;2]  [c(2-0)]2  ’ 


(25) 


This  means  that  a  nonspherical  body  for  general  values  of  aT  will  aquire  a  velocity  component  in  the  plane 
perpendicular  to  the  heat  current.  This  effect  is  apparent  in  the  final  behaviour  exposed  by  the  body  orbit  of 
the  coin  in  figure  1.  The  maximum  angle  between  the  heat  current  and  the  body  velocity,  i[>m,  occures  in  case 
of  purely  diffuse  reflexion  for  an  optimal  value  of  the  final  orientation  of  the  axis  of  symmetry  relative  to  the 
heat  current,  and  is  given  by 


ipM  =  arccos 


,\/T+ k 

'  2  +  k 


16.4°  (coin) 
12.0°  (needle) 


(26) 
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These  results  thus  predict  that  in  an  initially  localized  assembly  of  equatorially  symmetric  bodies  with  different 
orientations  relative  to  the  heat  current  will,  in  the  plane  perpendicular  to  q,  with  time  spread  without  bound. 
This  is  an  interesting  result,  but  perhaps  not  the  whole  story:  Going  to  higher  order  contributions  to  the 
torque,  a  term  of  the  order  of  Kn 2  will  recouple  the  angular  momentum  of  equatorially  symmetric  bodies  to 
the  temperature  gradient.  However,  since  in  the  present  context  Kn  -C  1,  this  coupling  is  clearly  weak. 


APPENDIX 


The  geometric  tensor  integrals  encountered  when  the  expressions  for  the  force  and  the  torque  are  integrated 
over  the  body  surface  are  listed  below: 

4|1,0)  =  ° 

7,(0,1)  =  53/2c(°.1)_/y, 


7S2°I  =  5  [K*  +  c^N<uNi2>\ 


K1’1)  =  s3/ 2 


*12 


41,%+41,1)A<iivi> 


m  =  s* 


Ann 


.(0.2), 


1  Sjih  +4  ’  ^  N<jt  Nj2  > 
(2,1)/?.  .  AT.nJ 2,!) 


I?:1].  =  S 3/2 


c^’S^Nj  +  cr>  (6idNi2  +  6i2jNh)  +  c^N^Ni 


j(2,2)  _  g2 

*i*2|iij2 

.(2,2) 


4  ’  +4  ^  )  +4  ’  ^ii»2-^jiAj2  +4 

+4  ’  *  (fihji ^i2Nj2  +  8i1j2Ni2Nj1  +  Si2j1Ni1  Nj2  +  Si2j2NitNjt) 

,(2 

'6 


+cg'2)NilNi2NhNh 


The  coefficients  in  the  expressions  above  are  given  by 


c(04)  =  j2 

g(2,0)  =  I  (3i7i  _ 

4X,1)  =  \J2 

4m)  =  I(3  j3-j2) 

4°’2)  =  i 

40,2)  =  |  (3j7  -  *75) 

=  l  (J2-J4) 

42,1)  =  5  (Js  -  *74) 

42.1)  =  I  (— 2v72  -  2J3  +  5  J4) 

42’2)  =  I  (-3 +  17Js  +  2^6  -  8J7  -  3Js  +  Jg) 
42’2)  =  |  (2*7s  -  2^6  -  3 Jr  -  3J8  +  ij9  +  J10) 
42’2)  =  I  (_ *75  +  2  J6  +  Jr  +  *7?  —  4Jg  +  Jx 0) 

42.2)  =  l  (-3*75  +  2Je  +  7  Jr  +  3Jg  -  4J9  -  5  Ji0) 
42’2)  =  I  (-3*75  +  2 J-6  +  3 J7  +  7Js  -  4J9  -  5  Ji0) 
42’2)  =  |  (-*7s  -  2  J6  -  Jr  -  *7g  -  8Jg  -  5Ji0) 

42,! 2)  =  |  (*75  +  2  J6  -  5 Jr  -  5Jg  -  5  J9  +  35Ji0) 
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In  these  coefficients  the  integrals  J\  —  J\9  are  given  by 


Intgrals 

Ji  =  S-1  Js  (N  •  n)2  dS 
J2  =  S~a/2fsx-NdS 
J3  =  5-3/2  fs  (x  •  n)  (n  •  N)  dS 
Ji  =  5-3/2  /  (x  •  N)  (n  •  N)2  dS 
=  S~2  Js  x.2dS 
J§  =  5“2  Jg  (x  •  n)2  dS 
J7  =  5-2  Js  (x  ■  N)2  dS 
Jg  =  S~ 2/  x2  (n  ■  N)2  dS 
J9  =  S~ 2  fs(x-  n)  (x  •  N)  (n  •  N)  dS 
Jio  =  S-2/s(x-N)2  (n  •  N)2  dS 

In  the  following  table,  the  values  of  the  coefficients  a±  —  a9  and  there  resp.  equatorial  symmetry  are  listed: 


Coefficient 

Value 

Eql. 

ai 

|  (4  -  3qt)  (3Ji  -  1) 

+ 

ai 

|  (8  +  7ro!T) 

+ 

0-3 

l  [8  -  (6  -  7T)  Or]  (3  Jx  -  1) 

+ 

(14 

i  {[8  +  (7T  -  2)  aT]  J2  [8  (6  7T)  aT\  J3) 

1  [(4  —  aT)  J2  —  (A  —  3ar)  J3] 

15  {I8  +  (tt  -  2)  aT]  J2-[8-(6- 1 r)  aT\  J73} 

- 

a5 

- 

OjQ 

- 

a7 

^  {[8  +  (7r  -  2)  aT]  Jb  -  [8  -  (6  -  7r)  aT]  J6} 

|  {[16  -  (10  -  2tt)  aT\  Jr,  -  [16  -  (12  -  2tt)  aT\  J6  -  [24  -  (12  -  3tt)  aT]  J7 
-  [24  -  (18  -  3tt)  aT\  J8  +  [48  -  (36  -  6tt)  ot]  J9} 

+ 

&8 

+ 

The  tensor  of  inertia,  Iij,  is  conveniently  expressed  according  to  (21)  with  the  choice 


I0  =  mBS, 


and  the  coefficients  61  and  b2  become 

bl  =  3m^S  Jv  62  =  2 ^8  jy  b2  -  3(N '  *)21 
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